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I. INTRODUCTION 

In this article, we report on the measurement of the 
branching fractions of the 22 decays of charged and neu- 
tral B mesons to D^D^K finaljtates (Table I): £>W 
is either a D°, D*°, D + or D* + , £)(*) is the charge con- 
jugate of £)(*) and K is either a K + or a K°. Both 
and £)M are fully reconstructed. Charge conjugate 
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reactions are assumed throughout this article. 

In the past, the values measured for hadronic decays 
of the B meson were in disagreement with the expec- 
tations based on the B semileptonic branching fraction 
due to the inconsistency originating from the number of 
charmed hadrons per B decay (charm counting) [1] . The 
b — >• ccs transition in B decays was believed to be dom- 
inated by B ->• D s A, B -> (cc) A, and B ^ S C X final 
states, where X represents any particles. However, it was 
realized [2] that an enhancement in the b — > ccs transi- 
tion was needed to resolve the theoretical discrepancy 
with the B semileptonic branching fraction. Buchalla 
et al. [2] predicted sizeable branching fractions for de- 
cays of the form B D^D^K (X). Experimental 
evidence in support of this picture soon appeared in the 
literature [3], including a study by BABAR using 76 fb _1 
of data where the Collaboration reported the observa- 
tions or the limits on the 22 decays B -> D^D^K [4]. 
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The aggregate branching fraction measurements were 
B(B° -> D^D^K) = (4.3 ± 0.3 ± 0.6)% and B(B+ 
£)(*)£)(*) iiT) = (3.5 ± 0.3 ± 0.5)%, where the first uncer- 
tainties are statistical and the second systematic. This 
result may be compared with the wrong-sign D produc- 
tion (b — > ccs transition containing a D meson) that 
BABAR studied using inclusive B decays to final states 
containing at least one charm particle [5]. The wrong- 
sign D production was found to be B(B° — > DX) = 
(10.4 ± 1.9)% and B(B~ — > DX) = (11.1 ± 0.9)%. In 
addition, BABAR found a value of the total charm yield 
per B decay consistent with the one derived from the 
semileptonic branching fraction, which solved the long- 
standing problem of the charm counting. 

Furthermore, D^D^K events are interesting for a 
variety of studies. These events can be used to inves- 
tigate isospin relations and to extract a measurement 
of the ratio of Y(4S) -> B + B~ and Y(4S) -> B°B° 
decays [6]. It was shown theoretically that the time- 
dependent rate for B° -> D^~D^+K° decays can be 
used to measure sin 2/3 and cos 2/3 [7]. BABAR used the 
mode B° -> D*-D*+K° s with 209 ftT 1 of data to per- 
form a time-dependent CP asymmetry measurement to 
determine the sign of cos 2/3, under some theoretical and 
resonant structure assumptions [8]. The Belle Collabo- 
ration also published a similar analysis [9] . Although the 
resonant states are not studied in our paper, it is worth 
recalling that many D^K and D^D^ resonant pro- 
cesses are at play in the studied decay channels. Using 
B -> D^D^K final states, BABAR and Belle observed 
and measured the properties of the resonances (2536), 
£> sJ (2700), V13770), and A(3872) [10-12]. 

The B -> D^D^K decays can proceed through ex- 
ternal IF-emission and internal IT-emission amplitudes, 
also called color-suppressed amplitudes. As Fig. 1 illus- 
trates, some decay modes proceed through only one of 
these amplitudes while others proceed through both. 

In this paper, we update with the full BABAR data sam- 
ple our previous measurement [4] of the branching frac- 
tions for the 22 B ->■ D^D^K and B -+ D^D^K+ 
decays. We benefit from several improvements with re- 
spect to this previous measurement: 

• the integrated luminosity used for this analysis is 
more than 5 times larger, 

• the track reconstruction and particle identification 
algorithms have been improved (in purity and effi- 
ciency), 

• the efficiency of the selection of signal events has 
been increased, 

• the fit uses a more accurate signal parametrization, 

• the peaking background is taken into account in the 
fit, 

• we use a method that is insensitive to the possible 
resonant structure in the final states. 



Measuring the 22 modes altogether allows to avoid biases 
in the branching fraction measurement by correctly tak- 
ing into account the cross-feed events, which are events 
from one mode being reconstructed as a candidate for 
another mode. 



TABLE I: The 22 B -> D w D w K decay modes. The 
modes B° ->■ D~D* + K° and B° ->■ D*~D + K° are com- 
bined together since they are not experimentally distinguish- 
able. The same applies to the modes B° -> D°D*°K° and 
B o ^ D*°D°K° which are also combined together. 



Neutral B mode 



Charged B mode 



B° ->■ D~D°K + 



B u 



D~D*°K + 
B° -> D*~D°K + 
B° -> D*~D* K + 
B° -> D-D+K 
B° -> D-D*+K° + D*-D+K° 

B° -> D*~D* + K° 

B° -> D°D°K° 

B° -> D°D*°K° + D*°D°K° 

B° -> D*°D*°K 



B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 
B+ 



D°D+K° 

D°D*+K° 

D*°D+K° 

D*°D*+K° 

D°D°K+ 

D°D*°K+ 

D*°D°K+ 

D*°D*°K + 

D-D+K+ 

D~D*+K + 

D*-D+K+ 

D *- D *+ K + 



II. THE BABAR DETECTOR AND DATA 
SAMPLE 

The data were recorded by the BABAR detector at the 
PEP-II asymmetric-energy e + e~ storage ring operating 
at the SLAC National Accelerator Laboratory. We an- 
alyze the complete BABAR data sample collected at the 
Y(AS) resonance corresponding to an integrated luminos- 
ity of 429 fb _1 , giving N B g = (470.9±0.1±2.8) x 10 6 BB 
pairs produced, where the first uncertainty is statistical 
and the second systematic. 

The BABAR detector is described in detail else- 
where [13]. Charged particles are detected and their mo- 
menta measured with a five-layer silicon vertex tracker 
and a 40-layer drift chamber in a 1.5 T axial magnetic 
field. Charged particle identification is based on the mea- 
surements of the energy loss in the tracking devices and 
of the Cherenkov radiation in the ring-imaging detector. 
The energies and locations of showers associated with 
photons are measured in the electromagnetic calorimeter. 
Muons are identified by the instrumented magnetic-flux 
return, which is located outside the magnet. 

We employ a Monte Carlo (MC) simulation to study 
the relevant backgrounds and estimate the selection effi- 
ciencies. We use EVTGEN [14] to model the kinematics 
of B mesons and JETSET [15] to model continuum pro- 
cesses, e+e~ — > qq (q = u,d,s,c). The BABAR detector 
and its response to particle interactions are modeled us- 
ing the GEANT4 [16] simulation package. 
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FIG. 1: Top left: external W-emission amplitude for the decays B + — > D^°D^ + K°. Top center: internal W-emission 
amplitude for the decays B + — > K + . Top right: external+internal W-emission amplitudes for the decays 

B + D(*)°D ( * )0 K+. Bottom row: same as top row respectively for B° -> D ( * )0 , B° -> D ( * )0 D M °K and 

B° -»• D^-D^+K . 



III. B CANDIDATE SELECTION 

We reconstruct the 5° and £> + mesons in the 22 
£)(*)£){*) X modes. The level of background widely varies 
among the signal channels, even within a specific B mode 
depending on the D meson decay type. A different opti- 
mization of the selection criteria is implemented for each 
of the final states. The optimization determines the selec- 
tion which maximizes S/y/S + B, where S and B are the 
expected number of events for the signal and for the back- 
ground in the signal region, based respectively on signal 
and background MC simulated events. The branching 
fractions for the computation of S are taken from our 
previous measurements of these channels [4] . 

We identify charged kaons using either loose or tight 
criteria depending on the decay mode. The loose crite- 
rion is typically 98% efficient with pion misidentification 
rates at the 15% level, while the tight criterion is 85% ef- 
ficient with a misidentification around 2%. We use only 
the meson when a neutral K meson is present in the 
final state. The candidates are reconstructed from 
two oppositely charged tracks assumed to be pions con- 



sistent with coming from a common vertex and having 
an invariant mass within ±9.5 MeV/c 2 of the nominal K® 
mass [17]. The displacement of the vertex in the 
plane transverse to the beam axis is required to be at 
least 0.2 cm. 

The 7T° candidates are reconstructed from pairs of pho- 
tons with energies E 1 > 30 MeV in the laboratory frame 
that have an invariant mass of 115 < m 77 < 150 MeV/c 2 . 

We reconstruct D mesons in the modes D° — > K~tt + , 
K-tt+tt , if-7r+7r-7r+, and /) -> K-tt+tt+. The K 
and 7r tracks are required to originate from a common 
vertex. The invariant masses of the D candidates are 
required to lie within ±2.5ctd of the measured D mass, 
where ctd is the D invariant mass resolution. This res- 
olution is measured to be 5.8 MeV/c 2 for D° — > K~tt + , 
9.5 MeV/c 2 for D° -> K-tt+tt , 4.7 MeV/c 2 for D° — > 
K-tt+tt-tt+, and 4.2 MeV/c 2 for D+ -> K~-k + -k+. To 
reduce the combinatorial background, for some of the B 
decays involving D° — > A"~7r + 7r°, we use the distribu- 
tion of events in the Dalitz plot of the squared invariant 
masses to 2 (A'~7t + ) x m 2 (A' _ 7r°) where we select events 
that are located in the enhanced regions dominated by 
the if* (892)+, A-*(892)°, and p(770)+ resonances [18]. 
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The D* candidates are reconstructed in the decay 
modes D*+ D°ir+, D*+ -> L>+7r°, D*° L>°7r 6 , 
and D*° — > D°7. The 7r° and the 7r+ candidates must 
have a momentum smaller than 450MeV/c in the Y(4S) 
rest frame, while the 7 energy in the laboratory frame 
must be larger than 100 MeV. The mass difference be- 
tween the D* + and D candidates is required to be within 
±3MeV/c 2 of the nominal value [17]. For D*° meson de- 
cays, the mass difference between the D*° and D° candi- 
dates is required to lie between 138 and 146 MeV/c 2 for 
D*° -> D° 7T° and between 130 and 150 MeV/c 2 for L>*° 
-> D° 7. 

The B candidates are reconstructed by combining a 
D^*\ a and a X candidate in one of the 22 modes. 
For modes involving two D° mesons, at least one of them 
is required to decay to K ~7r+, except for the decay modes 
D*-D*+K°, D* J D*+K+, and D*~D°K + , which have 
lower background and for which all combinations are ac- 
cepted. For modes containing a D* + meson, we look only 
to the decay D* + — >■ D°ir + , except for the modes contain- 
ing D*~ D* + , where we also reconstruct D* + -> D + ir°. A 
mass-constrained kinematic fit is applied to the interme- 
diate particles (D*°, D*+ , D°, D+, K°, tt°) to improve 
their momentum resolution. 

To suppress the continuum background, we remove 
events with R 2 > 0.3 (where R2 is the ratio of the sec- 
ond to zeroth Fox- Wolfram moments of the event [19]) 
and events with | cos(6>s)| > 0.9 (where 8b is the angle 
between the thrust axis of the candidate decay and the 
thrust axis of the rest of the event). 

Two kinematic variables are used to isolate the B- 
meson signal. The first variable is the beam-energy- 
substituted mass defined as 

f s/2+p .p B \ 2 m 

mEs = n — — J ~ bB| ' (1) 

where y/s is the e + e~ center-of-mass energy. For the mo- 
menta po, Pb and the energy E , the subscripts and B 
refer to the e + e~ system and the reconstructed B meson, 
respectively. The other variable is AE, the difference be- 
tween the reconstructed energy of the B candidate and 
the beam energy in the e + e~ center-of-mass frame. Sig- 
nal events have toes compatible with the known i?-meson 
mass [17] and AE compatible with MeV, within their 
respective resolutions. At this stage, we keep only events 
which satisfy toes > 5.20GcV/c 2 . 

We obtain a few signal B candidates per event on aver- 
age. When the final state contains no D* meson, we get 
1.0 to 1.3 candidates per event depending on the specific 
mode, 1.3 to 1.9 candidates per event for final states con- 
taining one D* meson, and 1.7 to 2.1 candidates per event 
when the final state contains two D* mesons (except for 
B+ -> D*°D*°K+ with 2.9 candidates per event). If 
more than one candidate is selected in an event, we retain 
the one with the smallest value of \AE\ ("best candidate 
selection"). According to MC studies, this criterion finds 
the correct candidate when this one is present in the can- 
didate list in more than 95% of the cases for final states 



with no D*° meson and more than 80% of the cases for 
modes with one or two neutral D* mesons. We keep only 
events with \AE\ < E c with E c varying from 7 MeV to 
56 MeV depending on the decay mode of the B and D 
mesons. The resolution on AE varies between 5.6 and 
14.3 MeV for modes with zero or one D*° meson in the 
final state, and between 11.6 and 19.5 MeV for modes 
containing two neutral D* mesons. 

The efficiency for signal events varies from 0.5% to 
22.2% depending on the final state (being typically in the 
5% — 10% range). The modes with the lowest efficiency 
are the ones containing one or two charged D* mesons. 

Figure 2 presents the AE and Toes distributions af- 
ter the complete selection is applied. The AE distribu- 
tions are presented for events in the signal region defined 
by Toes > 5.27GeV/c 2 and are shown without applying 
the best candidate selection. Signal events appear in the 
peak near AE ~ OMeV when reconstructed correctly, 
while the peak around —160 MeV is due to D*DK and 
DD*K decays reconstructed as DDK, and to D*D*K 
decays reconstructed as D*DK or DD*K. Both AE 
and toes distributions show a clear excess of events in 
the signal region. 



IV. FITS OF THE DATA DISTRIBUTIONS 

We present the fits used to extract the branching frac- 
tions. For each mode, we fit the toes distribution be- 
tween 5.22 and 5.30 GeV/c 2 to get the signal yield. The 
data samples corresponding to each B decay mode are 
disjoint and the fits are performed independently for each 
mode. According to their physical origin, four categories 
of events with differently shaped toes distributions are 
separately considered: D^D^K signal events, cross- 
feed events, combinatorial background events, and peak- 
ing background events. The total probability density 
function (PDF) is a sum of these contributions. Event 
yields are obtained from extended maximum likelihood 
unbinned fits. 



A. Signal contribution 

The shape of the signal is determined from fits to the 
Toes distributions of signal MC samples. A Crystal Ball 
function [20] (Gaussian modified to include a power-law 
tail on the low side of the peak) , Vs (toes ; tos , as , as , ns ) , 
is used to describe the signal (see Eq. (Al) in the Ap- 
pendix A). The parameters of this PDF are tos and as, 
the mean and the width of the Gaussian part, and as 
and ns, the parameters of the tail part. The signal yield, 
As, is determined from the fit to the data. 
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FIG. 2: Distributions of the AE variable (top plots) and of the mES variable (bottom plots) for the sum of all the B° — > 
D^D^K modes (left-hand plots) and for the sum of all the B + — > D^D^K modes (right-hand plots). The AE distributions 
are shown after the complete selection but before the choice of the best candidate and for 5.27 < mES < 5.29 GeV/ c 2 , and the 
mES distributions are shown after the complete selection, including the selection on the AE variable. 



B. Cross- feed contribution 

We call "cross feed" the events from all of the D^D^K 
modes, except the one we reconstruct, that pass the com- 
plete selection and that are reconstructed in the given 
mode. The cross-feed events are a non-negligible part of 
the tties peak in some of the modes, and the signal event 
yield must be corrected for these cross-feed events. 

We observe from the analysis of simulated samples that 
most of the cross feed originates from the combination 
of an unrelated soft ir° or 7 with the D° decayed from 
the D* + to form a wrong D*° candidate. The cross-feed 
proportion is often in the 10% range relative to the signal 
yield but can be comparable or larger than the signal 
contribution, especially for modes containing 
in the final state. To account for the cross-feed events, 
an iterative procedure, described in Sec. IV F, is used to 
extract the signal yields and the branching fractions. 

Cross-feed distributions for modes containing no 
D*° meson can be described by a Gaussian function 
•ppea ln s( mEg ; mcF, ccf) for the peaking part where mcF 
and acF are the mean and the width of the peaking com- 
ponent [Eq. (A2)]. For modes containing at least one neu- 
tral D* meson, the peaking component is described by 
a function V' C p ea lns (niFS', mcF, ccf> ^cf) which is able 
to model the tail at low mass [Eq. (A3)]. The parame- 
ters mcF and <Jqf represent the position of the maximum 
value and the width of the peak, and tcF represents the 



tail of the function. The nonpeaking part of the cross- 
feed contribution is described by an Argus function [21] 

7 , Qp npeakmg (mEs; mo, £cf) 5 where mo represents the kine- 
matic upper limit for the constrained mass and £cf is the 
Argus shape parameter [Eq. (A4)]. 
The total PDF for cross feed events is 

Vcf (m-ES ; m C F , c C f , ^cf , m , £cf ) = (2) 

^ypcaking ^ *p( f ) peaking . ^ynonpeaking ^ .pnonpeaking 
CF CF CF CF 

where V^) F poakmg represents either -pP° akm s or -p^peakmg 
depending on the number of neutral D* meson in the fi- 
nal state. The quantities A^ aking and N ™ npcakins are 
the numbers of events in the peaking PDF and in the 
nonpeaking PDF, respectively. The values of the param- 
eters of the cross-feed PDF are determined by fitting sig- 
nal MC mss distributions, except for the value of m 
which is fixed to 5.2892 GcV/c 2 . The cross-feed yield, 
iVcF = 7v£p npcakmg + iV£p akms , is also extracted from the 
fit. 



C. Combinatorial background contribution 

The combinatorial background events are composed of 
generic B decays and of continuum events, which ac- 
count, respectively, for about 88% and 12% of the to- 
tal number of background events. The combinatorial 
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background events are described by an Argus function 
T^cb^es; mo, £cb)j where £cb is the shape parameter 
[Eq. (A5)]. The parameter £cb is free to float in the fit 
to the data while m is fixed to 5.2892 GeV/c 2 . The yield 
for the combinatorial background, Acb, is also obtained 
from the data fit. 



D. Peaking background contribution 

We call "peaking background" the part of the background 
that is peaking in the signal region and that is not due to 
cross feed. To extract the peaking background, we fit the 
m ES distributions from generic MC samples e+e~ — > qq 
(q = u,d,s,c,b) satisfying the D^D^K selection and 
scale the results to the data luminosity. 

The simulated distribution is fitted with an Argus 
function describing the nonpeaking part and a Gaus- 
sian function T^peI^es; mpe, cpb) describing the peak- 
ing part, where mpe and ope are the mean and width of 
the Gaussian [Eq. (A6)]. The parameters mpe and ctpb 
are free to float in the fits to the simulated events, except 
for modes with nonconverging fits, where mpe is fixed to 
the B mass. These modes are B° -> D~ D+ K Q , B+ -> 
D°D+K°, B° -> D*-D°K+, B+ D*°_D + K° , B+ — > 
D°D*+K\ B° -> D*-D*°K+, B° -> D*° D*° K° , and 
B+ -> D*°D*°K+. The fit also returns the value of 
the peaking background yield, Apb, which is shown in 
Table II. Only the peaking part Vpb is used in the fit 
to the data, the nonpeaking part being included in the 
combinatorial background. 

E. Fits 

We fit the Toes distribution using the PDFs for the signal, 
for the cross feed, for the combinatorial background, and 
for the peaking background as detailed in the previous 
sections. The total PDF Vtot can be written as 

Vtot = A S X "Ps(wES;TOS,Cr S ,Q!S,«s) (3) 

+ A C f x Vcf (niES ; m C F , &cf , t cf , m , £ C f ) 
+ A C b x Pcb^es^o^cb) 
+ Apb x VpB(m ES ;m PB ,<T PB ). 

The free parameters of the fit are A s , m s , A CB , and 
£cb- All other parameters, except mo, are fixed to the 
values obtained from the simulation. For modes with 
low signal statistics in the data, namely B° — > D D°K° 7 
B" ->■ D°D* a K° + D*°D°K° and B a -> D*°D*°K°, we 
fix ms to the value obtained from the simulation. 

The free parameters are extracted by maximizing the 
unbinned extended likelihood 

p -N N n n 

£ =^r-n^ (4) 

where n is the number of events in the sample and A is 
the expectation value for the total number of events. 



F. Iterative procedure 

Because of the presence of cross-feed events, the fit for 
the branching fraction for one channel uses as inputs 
the branching fractions from other channels. Since these 
branching fractions are in principle not known, we em- 
ploy an iterative procedure. In practice, we perform the 
complete analysis for each B mode, using as a start- 
ing point the branching fractions measured by BABAR in 
Rcf. [4]. We obtain new measurements of the branching 
fractions that we use in the next step to fix the cross- 
feed proportion. We repeat this procedure until the dif- 
ferences between the actual branching fractions and the 
previous ones are smaller than 2% of the statistical un- 
certainty. Using this criterion, four iterations are needed. 
We keep the last iteration as the final result. 

G. Fit results 

The results of the fits are shown in Figs. 3 and 4, and are 
displayed in Table II. All the fits show a good description 
of the data. Although we perform an unbinned fit, we 
can compute a x 2 value using bins of 2.5MeV/c 2 width. 
We observe values of xV-^dof typically close to 1, with 
Adof = Abin — Afl oa t, where N a - m is the number of bins 
and Afl oa t is the number of floating parameters in the fit. 

V. BRANCHING FRACTION 
MEASUREMENTS 

A. Method 

In this paper, we measure the branching fractions of 
the 22 D^D^K modes, including nonresonant and 
resonant modes. It has been shown that D^D^K 
events contain resonant contributions. This was first 
reported by the BABAR Collaboration in Ref. [4] where 
it was observed that the three-body phase-space decay 
model does not give a satisfactory description of these 
decays. In a subsequent study [11], we showed the 
presence of .0^(2536), ^(3770) and A (3872) mesons in 
these final states. From Belle [10], we know that the 
D s j(2700) meson has a large contribution in the mode 
B + — > D°D°K + . This meson is expected to be present in 
DWdWK final states containing D°K+ and D+K°, as 
well as in the final states containing D*°K + and D* + K°, 
since it was recently seen decaying to D*K [22]. There is 
in addition the possibility of having unknown resonances 
in the D^D^K final states. Simulations of the known 
resonances indicate that the efficiencies for nonresonant 
modes and resonant modes are significantly different. 
This is due to the fact that the efficiency is not uniform 
across the phase space and that resonant events, depend- 
ing on the mass, the width and the spin of the resonance, 
populate differently the Dalitz plane. Ignoring this effect 
would introduce a bias of up to 9% in the total branching 
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FIG. 3: Fits of the mES data distributions for the neutral modes, B° — > D'*^D^K. The decay mode is indicated in the 
plots. Points with statistical errors are data events, the red dashed line represents the signal PDF, the blue long-dashed line 
represents the cross-feed event PDF, the blue dashed-dotted line represents the combinatorial background PDF, and the blue 
dotted line represents the peaking background PDF. The black solid line shows the total PDF. 



fraction for some decay modes In order to measure the 
branching fractions inclusively without any assumptions 
on the resonance structure of the signal, we estimate the 
efficiency as a function of location in the Dalitz plane of 
the squared invariant masses to 

(£)(*)£)(*)) X TO 
for the data. We use this efficiency at the event position 
in the Dalitz plane to reweight the signal contribution. 
To isolate the signal contribution event-per-event, we use 
the sPlot technique [23]. 

The sPlot technique exploits the result of the Toes 
fit (yield and covariance matrix) and the PDFs of this 
fit to compute an event-per-event weight for the signal 
category and background category. The PDF for the sig- 



nal category is Vs- For the background category, the 
three PDFs for the different components of the back- 
ground (cross feed, combinatorial events, and peaking 
background) are combined together to form one PDF: 



N C F X "PpF + ^CB X VcB + ^PB X V PB 



(5) 



where N B is the sum of the background yields, N B — 
N C f + N C b + Np B - The PDFs and the yields are the 
ones obtained from the results in Sec. IV. 

The sPlot weight for the signal category is defined as 



w s (i) 



Vs,sVs(i) + Vs,BVB(i) 
N s Vs(i) + N B V B (i) 



(6) 
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FIG. 4: Fits of the to-es data distributions for the charged modes, B + — ► D^*'D^*'K. The decay mode is indicated in the 
plots. Points with statistical errors are data events, the red dashed line represents the signal PDF, the blue long-dashed line 
represents the cross-feed event PDF, the blue dashed-dotted line represents the combinatorial background PDF, and the blue 
dotted line represents the peaking background PDF. The black solid line shows the total PDF. 



Here, i stands for the index of the event, with Vj(i) {j — 
S, B) corresponding to the value of the PDF for the event 
i. The quantity Vsj (j — S,B) is the covariance matrix 
element between the signal yield N$ and the yield Nj . 

The branching fraction for a specific D^*'D^*'K mode 
is given by 



B 



\ - w s (i) 
V N bb x (eB, 



sub/ 



where we define 



(eB, 



sub/ 



3 



x B. 



sub,j ■ 



(7) 



(8) 



The sum on j is over all the D subdecays of a particular 
£)(*)D^K mode. The term B su b,j is the product of the 
secondary branching fractions of the subdecay j: 



B. 



subj 



B 



DM 



x B DM x B 



A". 



(9) 



where B^ M ,B D (.), and Bk are the secondary branching 
fractions of the Z?M , £>W and K mesons [17] (with Bk 
1 for K + mesons). The quantity is the efficiency for 
the subdecay j at the Dalitz position of event i. In prac- 
tice, for a specific D^*>D^K mode with given D subde- 
cays {e.g. D° —> K + tt~ xfl"^ A'~7r + 7r°), the efficiency 
is obtained by using the specific simulated signal and di- 
viding the reconstructed signal by the generated signal 
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TABLE II: Number of events for the signal, N3, for the peaking background, Npb, and for the cross feed in the signal region, 
A^cp, and branching fractions in units of 10 -4 . The yields A^s and Ap B are defined on the whole m E s range, whereas A^p' is 
defined for tues > 5.27GeV/c 2 . The first uncertainties are statistical, and the second are systematic. The last column presents 
the significances including the systematic uncertainties. 

Mode As Apb Agp # Significance 

B° decays through external TV-emission amplitudes 
B° -»■ D-D {, K + 635 ± 47 99 ± 54 65 10.7 ± 0.7 ± 0.9 8.6a 



B° 


->■ D~D*°K + 




1116 ± 64 250 ± 69 137 34.6 ± 1.8 ± 3.7 


7.6a 


B" 


-¥ D*~D°K + 




1300 ± 54 93 ± 40 78 24.7 ± 1.0 ± 1.8 


12.6a 


7-l0 

B' 


-> D*~D*°K + 




1883 ± 63 31 ± 28 112 106.0 ± 3.3 ± 8.6 


11.4a 








decays through external±internal W-emission amplitudes 




7->0 

B' 


-¥ D~D + K" 




58 ± 10 8 ± 11 2 7.5 ± 1.2 ± 1.2 


5.1a 


dO 

B 


— > D D T A + D D^K 




422 ± 25 ± 12 7 64.1 ± 3.6 ± 3.9 


13.4a 


dO 

B 


— > D D T K 




511 ± 27 20 ± 13 5 82.6 ± 4.3 ± 6.7 


12.5a 








-B decays through internal TV-emission amplitudes 




B° 


-> D°D°K 




46 ±19 15 ± 19 19 2.7 ± 1.0 ± 0.5 


2.3a 


a 


L) JJ K + D UK 




126 ± 39 70 ± 39 147 10.8 ± 6.1 ± 3.0 


2.2a 


B 


. 7=7*0 r~i*0 r^0 

— > DDK 




170 ± 49 58 ± 31 231 24.0 ± 5.5 ± 6.7 


2.2a 








B + decays through external TV-emission amplitudes 




B+ 






237 ± 30 40 ± 23 16 15.5 ± 1.7 ± 1.3 


6.6a 


B T 


— )■ L> iJ 




233 ±19 9 ± 10 17 38.1 ± 3.1 ± 2.3 


10.7a 


B+ 


-> D*°D + K° 




164 ± 37 48 ± 33 95 20.6 ± 3.8 ± 3.0 


3.3a 


B+ 


-> D*°D* + K° 




308 ±28 11 ± 12 113 91.7 ± 8.3 ± 9.0 


7.5a 








decays through external±internal TV-emission amplitudes 




B+ 


-> _D° J D°A'+ 




901 ± 54 173 ± 77 153 13.1 ± 0.7 ± 1.2 


8.6a 


B+ 


-»• D°D*°K+ 




2180 ± 74 92 ± 50 409 63.2 ± 1.9 ± 4.5 


12.5a 


B+ 


-> D*°D°K+ 




745 ± 60 61 ± 26 724 22.6 ± 1.6 ± 1.7 


8.3a 


B+ 


-> D*°D*°K + 




3530 ± 141 186 ± 65 928 112.3 ± 3.6 ± 12.6 


6.8a 








B + decays through internal TV-emission amplitudes 




B+ 


-> D-D+K+ 




60 ± 15 35 ± 20 7 2.2 ± 0.5 ± 0.5 


2.8a 


B+ 


-> D-D*+K+ 




91 ± 13 2 ± 7 10 6.3 ± 0.9 ± 0.6 


6.7a 


B+ 


-> D*-D+K+ 




75 ± 13 15 ± 9 6 6.0 ± 1.0 ± 0.8 


5.1a 


B+ 


-> D*~D* + K + 




232 ± 23 30 ± 14 31 13.2 ± 1.3 ± 1.2 


7.4a 



in the Dalitz plane m 2 (D^D^) x m 2 (D<-^K), which is 
divided in 15 x 15 bins for the operation. The size of the 
bins is roughly 0.52 x 0.77GcV 2 /c 4 , 0.46 x 0.68GcV 2 /c 4 , 
and 0.38 x 0.58GcV 2 /c 4 for decay modes with no D* , 
one D* and two D* mesons respectively, depending on 
the available phase space. Neighboring bins are added 
together if one bin contains fewer than 10 events in the 
reconstructed Dalitz plane. The signal is simulated as- 
suming a flat (phase space) distribution in this Dalitz 
plane. 

The statistical uncertainty on the branching fraction 
is given by [23] 



w s (i) 



n bb x (£#sub); 



(10) 



B. Validation 

The analysis is validated at all stages by use of MC sam- 
ples. These samples consist of a mixture of continuum 
events and generic B decays containing the D^D^K 
signals with branching fractions close to the ones mea- 
sured in our previous result [4] . As a preliminary remark, 
it has to be noted that the analysis technique, including 
the selection optimization and the procedure for the fit 
and for the branching fraction measurement, is first de- 
termined solely on MC simulations ( "blind" analysis) . 

First, we show that the fit is able to find the true num- 
ber of simulated signal events within a la interval for the 
22 modes, where a is the statistical uncertainty reported 
by the fit. 

Furthermore, the sPlot method is tested on simulated 
samples. It is shown that this technique is able to tag 
true MC signal events with very good performance. A 
feature of the sPlot method is that the sum of sPlot 
weights for a given category is equal to the yield of this 
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category [23]. We determine that the sum of sPlot signal 
weights for the MC signal events is equal to the number 
of simulated MC signal events with a relative difference 
smaller than 1.5% for the majority of the modes. We also 
check that the sum of sPlot signal weights for the MC 
background events is compatible with zero as expected. 

Finally, we perform the measurement on MC simula- 
tions and find that the analysis is able to find the branch- 
ing fractions set in the simulation within a la interval for 
most of the 22 modes, where a is the total uncertainty 
on the branching fraction (combining in quadrature sta- 
tistical and systematic uncertainties). We also test the 
iterative procedure by randomizing the initial branching 
fractions and check that the branching fractions are con- 
verging to the expected values after a few iterations. 

C. Measurement 

For each event, we obtain the sPlot weight as well as 
the efficiency at its Dalitz position. Using Eq. (7), we 
compute the branching fraction for each of the B modes. 
We present these results in Table II. We assume equal 
B°B° and B + B~ production [17]. 

VI. SYSTEMATIC UNCERTAINTIES 

We consider several sources of systematic uncertainties 
on the branching fraction measurements. Their contri- 
butions are summarized in Table III. 

(a) We fix the width of the signal PDF from the value 
obtained in the fit to the signal MC sample. To esti- 
mate the systematic uncertainties originating from this 
choice, we repeat the fit with the width free to float for 
modes with high significance, namely B° — > D~ D* + K°+ 
D*-D+K\ B° -> D*-D*+K°, B° -> D*-D*°K+ and 
B+ -> D°D*+K°. The difference observed between 
the width from the data and from the MC events is 
roughly equal to 0.1MeV/c 2 . Using this number, we re- 
peat the fits for all modes adding ±0.1MeV/c 2 to the 
value of the width. The difference with the nominal 
branching fraction gives the systematic contribution as- 
sociated to the signal shape. In addition, as mentioned in 
Sec. IV E, three D^D^K modes with low signal statis- 
tics have their PDF mean fixed to the value obtained for 
the simulation. We repeat the fit using the PDF mean 
obtained from a mode with a large statistics (namely 
B° -> D*~ D* + K°) and take the difference in branching 
fraction as the systematic uncertainty. These two con- 
tributions of the systematic uncertainties are added in 
quadrature. 

(b) The cross-feed determination introduces systematic 
uncertainties of which two sources are identified. First, 
we use an alternate function for the cross-feed PDF [us- 
ing a nonparametric function rather than Eq. (2)], which 
gives relative systematic uncertainties below 1%. Second, 
the cross-feed branching fractions and their uncertainties 



are known from the results of this analysis. To estimate 
the systematic uncertainties coming from this effect, we 
repeat the measurement applying ±lcr of the statistical 
uncertainty on each cross-feed contribution to a given 
mode. These different contributions for each cross-feed 
mode are then combined quadratically. 

(c) The peaking background contributions are fixed 
from fits to the background MC simulation, using a Gaus- 
sian PDF. In these fits, the three parameters, namely, 
the number of events, the mean, and the width, are cor- 
related. We generate several sets of these parameters 
based on the covariance matrix of the fits and recom- 
pute the branching fractions for each of these sets. From 
the distribution of the branching fractions, we extract 
the systematic uncertainties originating from the peak- 
ing background. Another systematic effect arises from 
the fact that we use the MC after having scaled it to the 
data luminosity (using the total number of MC events 
passing the selection, the number of BB pairs, and the 
cross-section of e+e~ — > qq, where q — u,d,s,c). We 
estimate the data-MC agreement by computing the ra- 
tio of number of events for data and simulation for toes 
between 5.22 and 5.25 GcV/c 2 . We rescale the peaking 
background events using the ratio found in the specific 
mode (0.9 in average) and repeat the branching fraction 
measurement. The difference with the nominal branch- 
ing fraction gives the systematic uncertainty related to 
this effect. We combine these two sources of uncertain- 
ties in quadrature. Given the difficulty of estimating the 
peaking background, this is the dominant systematic un- 
certainty for most of the D^D^K modes. 

(d) The systematic uncertainty associated with the as- 
sumption of a fixed value of the end point in the Argus 
function is estimated by repeating the fit and letting the 
end point free to vary in the physical region between 
5.288 and 5.292 GeV/c 2 to account for possible varia- 
tions in the beam energy measurement. The difference 
in the branching fraction between this fit and the nom- 
inal fit gives the systematic contribution related to the 
combinatorial background. 

(e) We investigate the fit procedure performing a large 
number of test fits to MC samples obtained from the 
PDFs fitted to the data and look for the presence of pos- 
sible bias in the number of signal events. We observe 
that the biases are in most cases smaller than 10% of the 
statistical uncertainty. We do not correct these small bi- 
ases but take them into account in the total systematic 
uncertainties. 

(f ) An iterative procedure is performed to compute the 
branching fractions. We check this procedure on the MC 
simulation, where the results should not depend on the 
procedure used. The difference between the iterative and 
noniterative methods is small but non-negligible in some 
cases. We take the relative difference as the systematic 
contribution on the data due to the iterative procedure. 

(g) The limited Monte Carlo statistics induce an un- 
certainty on the computation of the signal efficiency. We 
use an efficiency mapping in the Dalitz plane with 15 x 15 
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TABLE III: Summary of the absolute systematic uncertainties on the branching fractions for each D^D^K mode (in units of 
1CP 4 ). The values listed in this table correspond to the systematic uncertainties associated with the signal shape (a), the cross- 
feed contribution (b), the peaking background (c), the combinatorial background (d), the fit bias (e), the iterative procedure 
(f), the limited MC statistics (g), the number of bins of the Dalitz plane (h), the particle detection efficiency (i) and finally 
the secondary branching fractions and the number of BB pairs (j). The letters in parenthesis refer to the specific paragraph in 
Sec. VI. The last column presents the total systematic uncertainties. 



Mode 


Signal Cross Peaking Comb. 


Fit 


Iter. MC 


Bins 


Particle 


BF + 


Total 






shape feed back. back, 
(a) (b) (c) (d) 


bias 
(e) 


proc. stat. 
(f) (g) 


(h) 


detection 

(il 
W 


AT — 

BB 

fil 
\J) 


syst. 






B° decays through external W 7 - 


emission 


amplitudes 










B' 


->■ D D°K + 


0.2 0.1 0.5 0.0 


0.0 


0.0 0.2 


O.z 


U.o 


0.3 


U.9 


dO 

B 


— > D~ D*°K + 


0.7 0.1 2.4 0.0 


0.0 


0.1 1.3 


1.1 


1.9 


1.1 


3.7 


dO 

B 


— > D D 


0.5 0.0 0.7 0.0 


0.0 


0.0 0.5 


0.5 


1.3 


0.8 


1.8 


D 


—^U JJ K 


2.0 0.2 1.6 0.1 


0.0 


0.2 1.8 


3.8 


6.4 


2.9 


8.6 






B decays through external+internal VF-emission amplitudes 










DO 
B 


— > JJ JJ K 


0.1 0.0 1.0 0.0 


0.2 


0.0 0.1 


1 


9 


n a 


1 9 


R° 


—t JJ JJ i\ -+- JJ JJ K 


0.9 0.0 1.3 0.1 


0.1 


0.0 1.3 


0.9 


2.4 


2.2 


3.9 


R° 


— r JJ JJ J\ 


1.2 0.0 2.0 0.0 


0.0 


0.0 3.2 


2.1 


4.2 


2.9 


6.7 






B° decays through internal W- 


emission 


amplitudes 










B' 


— > D'D K' 


0.2 0.0 0.4 0.1 


0.1 


0.0 0.1 


0.2 


U. 1 


U.l 


u.o 


R° 




0.3 0.4 3.4 0.1 


0.5 


0.1 0.5 


0.3 


0.5 


0.3 


3.6 


J j 


— T i t i t l\ 


0.9 1.6 5.2 0.2 


0.8 


1.8 1.1 


2.6 


1.7 


0.6 


6.7 






B + decays through external 


-emission 


amplitudes 










B+ 


-»• D°D+K° 


0.4 0.0 0.9 0.1 


0.0 


0.0 0.3 


0.2 


0.5 


0.5 


1.3 


B+ 


-> D°D*+K° 


0.5 0.3 1.0 0.1 


0.1 


0.0 0.8 


0.3 


1.4 


1.0 


2.3 


B+ 


_j. d*°D+K° 


0.5 0.3 2.5 0.3 


0.3 


0.3 0.8 


0.5 


1.0 


0.7 


3.0 


B+ 


-> D*°D* + K° 


1.6 1.3 2.5 0.1 


0.1 


0.4 3.1 


5.5 


4.9 


2.4 


9.0 






B + decays through external+internal W-emission amplitudes 








B+ 


-> _D° J D°K+ 


0.3 0.0 1.0 0.0 


0.0 


0.0 0.2 


0.2 


0.6 


0.3 


1.2 


B+ 


-> D°D*°K+ 


1.2 0.2 1.0 0.0 


0.1 


0.3 0.9 


1.4 


3.6 


1.6 


4.5 


B+ 


-> D*°D°K+ 


0.4 0.2 0.5 0.0 


0.1 


0.3 0.6 


0.5 


1.3 


0.6 


1.7 


B+ 


-> D*°D*°K+ 


2.0 0.6 2.8 0.7 


0.2 


1.3 4.2 


6.2 


9.0 


3.0 


12.6 






B + decays through internal W- 


■emission 


amplitudes 










B+ 


-> D~D+K+ 


0.1 0.0 0.5 0.0 


0.1 


0.0 0.0 


0.0 


0.1 


0.1 


0.5 


B+ 


-»■ D-D*+K+ 


0.1 0.0 0.4 0.0 


0.0 


0.0 0.2 


0.3 


0.3 


0.2 


0.6 



B + ^D*-D + K + 0.1 0.0 0.6 0.1 0.0 0.0 0.3 0.1 0.3 0.2 0.8 

B+ -> D*-D*+K + 0.2 0.0 0.5 0.0 0.1 0.0 0.6 0.2 0.7 0.4 1.2 



bins. To take into account this uncertainty, we generate 
several efficiency mappings, where in each bin we vary 
the nominal efficiency according to the efficiency uncer- 
tainty distribution. We obtain a distribution of branch- 
ing fractions that we employ to determine the systematic 
contribution. 

(h) We extract the efficiency in the Dalitz plane from a 
15x15 bin mapping. We vary the numbers of bins from 
lxl bin to 20 x 20 bins and recompute the branching frac- 
tions. This test is performed on MC simulation contain- 
ing signal and background events, where the D^D^K 
signal is purely nonresonant. In this case, the results 
should not depend on the number of bins since no res- 
onant states are present. The maximal difference with 
the nominal branching fraction is taken as the system- 
atic uncertainty 

(i) From the differences in the reconstruction and parti- 



cle identification efficiencies for the data and MC control 
samples, we derive systematic uncertainties of 0.2% per 
charged track, 1.7% per soft pion from D* decays, 1.2% 
per Kg, 3% per 7r°, and 1.8% per single photon. Addi- 
tionally, the systematic uncertainties for the K + identi- 
fication are ranging from 2% to 4% (in total) depending 
on the mode. 

(j) Finally, the uncertainties on the £)W and K® 
branching fractions [17] are accounted for. The total sys- 
tematic uncertainty also takes into account the number 
of B mesons in the data sample, which is known with a 
0.6% uncertainty. 

Table III shows a summary of the systematic uncer- 
tainties. The uncertainties from the different contribu- 
tions are added together in quadrature to give the total 
systematic uncertainty for a specific mode. 
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VII. RESULTS 

The final results on the data using the full BaBAR data 
sample can be found in Table II. In this Table, the 
quantity N^,p is the number of cross-feed events in the 
signal region (i.e. integrating the cross-feed PDF for 
toes > 5.27GeV/c 2 ) determined from the MC simula- 
tion scaled to the data luminosity using the branching 
fractions measured in this analysis (this includes peak- 
ing and nonpcaking cross- feed contributions). We in- 
dicate the significances (including systematic uncertain- 
ties) of the observations. To compute these significances, 
we repeat the fits using no contribution from the sig- 
nal. We compute the statistical significance 5 s tat calcu- 
lating PROB[21n(£ signa i/£o), Ndof], where £ s i gna i (Co) is 
the maximum of the likelihood with (without) the signal 
contribution, A^of is the number of free parameters in the 
signal PDF (two here) , and PROB is the upper tail proba- 
bility of a chi-squared distribution, converting this proba- 
bility into a number of standard deviations. We then take 
the systematic uncertainty into account by smearing S^tat 
by use of a Gaussian with a width equal to the system- 
atic uncertainty: 5 sta t+syst = S sta .t/^J (1 + er s 2 yst /(7 2 tat ), 
where <7 s t a t and a sys t are, respectively, the statistical and 
systematic uncertainties on the branching fraction mea- 
surement. 

We check isospin invariance using the D^D^K de- 
cays. Assuming isospin invariance in the B decay, inter- 
changing the u and d quarks in the Feynman diagrams of 
Fig. 1 should not modify the amplitude values. Table IV 
presents the ratios of the modes which are related by 
isospin symmetry. In the ratio of the branching fractions, 
all factors cancel except the amplitudes and the B° /B + 
lifetimes (neglecting the small mass differences between 
neutral and charged states for the B, D* , D, and K 
mesons) . We multiply the ratios of the neutral to charged 
branching fractions, r, by the ratio of the charged to neu- 
tral B meson lifetimes, t b +/t b o = 1.071 ± 0.009 [17]. 
The uncertainty on these values reported in Table IV 
combines the statistical and systematic uncertainties of 
our measurement, as well as the uncertainty on the life- 
time ratio. The values of r x t b + /t b o should be equal 
to unity if isospin invariance is verified. Although some 
values are compatible with this equality, for some others 
we observe discrepancies up to 2.6a (where a is the 68% 
standard deviation). This result is obtained assuming 
equal production of B° and B + mesons. 



VIII. CONCLUSION 

We have analyzed 471 x 10 6 pairs of B mesons produced 
in the BABAR experiment, and studied the exclusive de- 
cays of B°/B°, B ± to and B°/B°, B ± to 
D(*)D^K°. We measure the branching fractions for 
the 22 modes (see Table II). Some of the modes have 
been observed for the first time here: B + — > D*°D°K + 



(8.3cr), B+ -> D*~D*+K+ (7.4(7), B+ -> D-D*+K+ 
(6.7cr), B+ — > D a D+K a (6.6cr), and B Q ->■ D~D + K° 
(5.1a). In addition, we show evidence for the mode 
B+ -> D* Q D+K° (3.3(7) for the first time. We also report 
the observation of some of the color-suppressed modes, 
namely B+ -> D*-D*+K+ (7Aa), B+ -> D-D*+K+ 
(6.7cr), and B+ -> D*~ D+ K+ (5.1a). The other color- 
suppressed modes are seen with a lower significance: 
B+ -> D-D+K+ (2.8a), B° -> D a D a K a (2.3a), B° -> 
D*°D*°K° (2.2a), and B° -> D°D*°K° + D*°D°K° 
(2.2a). 

Summing the 10 neutral modes and the 12 charged 
modes, we measure that D^D^K events represent 
(3.68 ±0.10 ±0.24)% of the B° decays and (4.05 ± 0.11 ± 
0.28)% of the B + decays, where the first uncertainties 
are statistical and the second systematic, taking into 
account the correlations amongst the systematic uncer- 
tainties. These decays do not saturate the wrong-sign 
D production and account roughly for one third of this 
production. This result implies that probably decays 
of the type B -> D^D^K(n-ir) (with n > 1) have 
a non-negligible contribution to the b — > ccs transition 
(for example through the decays B -» D^D^K* or 
B — > D^D**K where D** is an excited D meson other 
than D*+ and D*°). 

The results obtained here are found to be in satis- 
factory agreement with those of the previous study by 
BABAR using 76 fb _1 [4] and supersede these previous 
measurements. Our branching fraction measurement of 
the mode B° — > D*~ D* + K° is found in good agreement 
with the values reported in Ref. [8] and Ref. [9], and su- 
persedes our previous result [8]. However, our branching 
fraction measurement of the mode B + — > D°D°K + is in 
disagreement at a 2.1c level with the Belle result [10]. 

We believe that the discrepancy with Ref. [10] and the 
fact that the branching fractions measured here are al- 
most systematically lower than the ones in Ref. [4] (al- 
though most of the time compatible) are due to the fact 
that in the present work we employ a more accurate 
parametrization of the signal toes distribution in the fit 
and that we take into account both the cross-feed and the 
peaking background contributions. In Ref. [4], only cross 
feed was accounted for but the peaking background was 
not considered due to the lower statistics. In addition, 
the efficiency correction used in obtaining the branching 
fractions accounts for the presence of resonant interme- 
diate states in the data. 

Finally, from neutral to charged B meson ratios of the 
branching fractions, assuming equal B°B° and B + B~~ 
production and taking into account the B meson life- 
times, we note that some mode ratios respect the isospin 
invariance, while some others show discrepancies up to 
2.6cr. 
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TABLE IV: Ratios of neutral to charged branching fractions. The second column shows the ratio r, where the first uncertainty 
is statistical and the second is systematic. The third column shows this value multiplied by the ratio of the charged to neutral 
B meson lifetimes, where the error includes all uncertainties. 



Mode 


T 


7" X TnX / Ton 


B(B° -»• D-D°K + )/B{B+ -> D°D + K°) 


0.69 ± 0.09 ± 0.08 


0.74 ± 0.13 


B(B° ->■ D-D M, K+)/B(B+ -> D°D*+K°) 


0.91 ± 0.09 ± 0.11 


0.97 ± 0.15 


B(B° -»• D*-D i) K + )/B{B+ -> D* i] D + K°) 


1.20 ± 0.23 ± 0.20 


1.28 ± 0.32 


B(B° ->■ D*-D*°K+)/B{B+ -> D*°D*+K°) 


1.16 ± 0.11 ± 0.15 


1.24 ± 0.20 


B(B° -»■ D-D+K°)/B(B + -> D°D°K+) 


0.57 ± 0.10 ± 0.10 


0.61 ± 0.15 


B(B° -> D~D*+K + D—D+K ) 
B(B+ -»■ D°D*°K+) +B(B+ -> D*°D°K+) 


0.75 ± 0.05 ± 0.07 


0.80 ± 0.09 


B(B° -> D—D* + K Q )/B{B+ ->• D*°D*°K+) 


0.74 ± 0.05 ± 0.10 


0.79 ± 0.12 


B(B° -»• D D {, K°)/B{B + -> D~D+K+) 


1.20 ± 0.53 ± 0.36 


1.28 ± 0.69 


B(B° -»■ D°D r0 K° + D*°D°K°) 


0.88 ± 0.27 ± 0. 31 


0.94 ± 0.44 


B(B+ -»■ £>-£>*+#+) + 23(B+ ->■ 7> -£>+#+) 


-»• D*°D*°K°)/B(B + -> D*-D*+K+) 


1.81 ± 0.45 ± 0.53 


1.94 ± 0.75 
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For modes containing at least one neutral D* meson, 
the peaking component is represented by an empirical 
function describing an asymmetric peak: 



p^pea ln S( mEg; mcF; cr CF , i CF ) 



(A3) 



exp 



ln 2 (l + MM 



TO C F 



2^CF 



0"CF 



t 2 



Appendix A: Fit PDF expressions 

We give the expressions of the PDFs introduced in 
Sec. IV (along with their parameters) and used to fit 
the toes distribution in the data. 



1. Signal PDF 

The signal PDF is given by 

Vs (toes ; m s , o-s , «s , n s ) 



cxp[-(to es - to s ) 2 /(2o|)] 

(ff)"s e xp(-al/2) 
^ ((ms - to E s)/cts + - a s ) r 



(Al) 

toes > to s - as os 
toes < to s - aso-s. 



In this equation and in the following, we omit the factor 
that normalizes the PDF to unity. 



2. Cross- feed PDF 

The cross-feed PDF for modes containing no D*° meson 
is described for the peaking part by 



ppea m S( mEs;mcF)( j CF ) = 

cxp[-(m ES - to C f) 2 /( 2<t cf)]- 



(A2) 



where A = sinh(i CF Vm 4) / (t CF vln 4) . This function ap- 
proaches a Gaussian function when the parameter i CF 
vanishes. 

The nonpeaking part PDF is 



m ESy /l - (toes/too) 2 x cxp[-^cf(1 - (to E s/to ) 2 )]. 
3. Combinatorial background PDF 

The combinatorial background PDF can be expressed as 



"Pcb (toes; mo, Ccb) = toes V 1 - (to E s/to ) 2 (A5) 
xcxp[-£ C b(1 - (to.es/toq) 2 )]. 



4. Peaking background PDF 

The peaking background PDF is given by 

"Ppb(to E s;topb,o-pb) = (A6) 
cxp[-(to es - mp B ) 2 /(2o-p B )]. 



